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Crystal point group symmetry is shown to protect Majorana fermions (MFs) in spinfull supercon- 
ductors (SCs). We elucidate the condition necessary to obtain MFs protected by the point group 
symmetry. We reveal that superconductivity in Sr2Ru04 hosts a topological phase transition to a 
topological crystalline SC, which accompanies a d- vector rotation under a magnetic field along the 
c-axis. Taking all three bands and spin-orbit interactions into account, symmetry-protected MFs 
in the topological crystalline SC are identified. Detection of such MFs provides evidence of the 
d-vector rotation in Sr2Ru04 expected from Knight shift measurements but not yet verified. 

PACS numbers: 



There has been a recent interest in the realization of 
Majorana fermions (MFs) in topological superconductors 
(SCs) [lHil. While it is known that spin-triplet SCs 
can host topological superconductivity [4 8], a recent 
breakthrough indicates that conventional ,s-wave super- 
conducting states may also support the top olog ical phase 
in the presence of spin-orbit interactions [9|— 1_1 



131. The fea- 



tures of s-wavc superconductivity and its possible appli- 
cation to fault-tolerant topological quantum computation 
have stimulated both theoretical and experimental activ- 
ities. The s-wave superconducting scheme has been ap- 
plied to a wide class of condensed matter systems from 
supcrfluids_of cold atoms to superconducting nanowire 
-H. 



14 



systems 

While these developments are based on topological 
classifications using the general symmetries of time- 
reversal and charge conjugation 2l|, systems often have 
other symmetries specific to their structures such as 
translational, rotational, and point group symmetries. 
Interestingly, additional symmetries can give rise to a 
nontrivial topology of the bulk wave functions and gap- 
less states on the boundaries [22rl26j . Although these spe- 
cific symmetries are microscopically sensitive to a small 
disturbance, recent studies of topological crystalline insu- 
lators have shown that if the symmetries are preserved on 
average, then the existence of gapless boundary states is 
rather robust (27rl29| . Therefore, it is expected that the 
symmetry-protected topological phase can provide an al- 
ternative platform for realizing MFs. 

In this Letter, we clarify how the crystal point group 
symmetry may protect the existence of MFs in SCs. We 
focus here on spinfull SCs with point group symmetry. 
Our arguments are directly applicable to many uncon- 
ventional SCs, most of which are spinfull. 

Using a representation of the gap function for the point 
group symmetry, we first elucidate the condition neces- 
sary to obtain MFs protected by the point group symme- 
try. We demonstrate that when this condition is satis- 
fied, spinfull SCs can be separated into a pair of spinless 
topological SCs, each of which supports MFs. If the con- 



ditions are not satisfied, however, the system reduces to 
a pair of states that are topologically of the same class as 
quantum Hall states, and thus they only support Dirac 
fermions at most. 

As a concrete example, we apply these arguments to 
the two-dimensional SC in Si^RuC^. It is argued that 
Si"2Ru04 hosts a profound topological phase transition 
from a chiral topological SC to a topological crystalline 
SC accompanied by d-vector rotation under a magnetic 
field parallel to the c-axis. We identify the symmetry- 
protected MFs by taking all three bands and spin-orbit 
interactions into account. Detection of the symmetry- 
protected MFs provides a distinct signature of the d- 
vector rotation in Sr2Ru04. Such a rotation is expected 
from Knight shift measurements but is yet to be verified. 

We begin with a description of two-dimensional SC 
based on the Bogoliubov de Genne (BdG) Hamiltonian 

« = E fc *l«(fe)*fc/2: 

n ^ - ( ) •** = (1) 

where Cksi is the annihilation operator of electrons with 
momentum k = (k x , k v ) and spin s, I denotes the orbital 
degrees of freedom of the electron, £(k) is the Hamil- 
tonian of the normal state, and A(fc) is the gap func- 
tion of the SC. Here, the spin (s =t,-l) and orbital 
(I = 1, ■ • • ,7V) indices are implicit in £(k) and A(fe). 
Both £(k) and A(fc) are 27V x 27V matrices. Assuming 
that the normal state has a mirror symmetry with re- 
spect to the xy-plane as M xy £{k)M xy = £ (fe), where 
A4 xy is the unitary matrix of the mirror reflection, we 
demonstrate how this mirror symmetry ensures topolog- 
ically stable MFs. 

We note that the superconducting state retains the 
mirror symmetry if the gap function A(fe) is even or odd 
under the mirror reflection, M xy A(k)M xy = ±A(fc). 
In the former (latter) case, 'H{k) is invariant under the 
mirror reflection 



M xy H{k)M xy = H{k), 



(2) 



2 



with M xy = M+ y (M xy ) given by, 



M 



M X y 

±M xy 



(3) 



In both these cases, H(k) commutes with either of M^ y 
(which we refer to simply as Ai xy in the following). 
Therefore, TL(k) is block diagonal in the diagonal ba- 
sis of M. xy . Each block-diagonal subsector has a definite 
eigenvalue of the matrix M. xy . 

The mirror Chern number ^(A) is defined as a Chern 
number of the subsector; using the negative energy states 
\u^(k)) in the subsector with the eigenvalue A of M xy , 
the gauge field in the momentum space is introduced as 



Am = z£ 



£„<(r 



i*(k)\dk a u*(k)). The mirror Chern 



number is then given by 



1 f* r n 
v{X) = —J J dk x dk y F x , 



(4) 



where J- x is the field strength of A x , and the integration 



is performed over the first Brillouin zone, —n < k XiV 



< 7T. 



The mirror Chern number introduced here is a natu- 
ral generalization of that used to characterize topological 
crystalline insulators [22| . Therefore, as well as topolog- 
ical crystalline insulators, the mirror Chern number en- 
sures the existence of gaplcss boundary states as long as 
the mirror symmetry is preserved macroscopically. How- 
ever, as we show, there is an important difference: to 
stabilize the MFs using the mirror symmetry, an addi- 
tional requirement for M xy is needed. 

This additional requirement originates from the sym- 
metry specific to SCs. Because of the self-conjugate prop- 
erty of the quasiparticle field ^fc in Eq. (JXJ) , 



C = 



1 

1 



(5) 



the BdG Hamiltonian has the special symmetry 
CU{k)tf = -H*(-k); i.e., particle-hole symmetry. The 
self-conjugate property of the particle-hole symmetry is 
the origin of the Majorana nature of topologically pro- 
tected gapless states. Therefore, to obtain MFs in a sub- 
sector of the system, particle-hole symmetry, which is 
closed in the subsector, is essential. A subsector of H(k) 
with a definite eigenvalue of A4 xy , however, does not al- 
ways have particle-hole symmetry within the subsector 
because this symmetry can exchange a pair of subsectors 
with different eigenvalues. Thus, only subsectors with 
particle-hole symmetry can support MFs that are pro- 
tected by the mirror symmetry. 

We now elucidate the condition for a subsector with an 
eigenvalue A of M xy to host its own particle-hole sym- 
metry. Since the particle-hole symmetry maps a state 
\u(k)) with the eigenvalue A to C\u*(—k)), the condition 
is derived so that the mapped state C\u*(—k)) has the 



same eigenvalue A as the original one. We find that this 
leads to 



CM XV & = X 2 Mt 



(6) 



In other words, only for M. xy satisfying Eq. ([6]) can one 
can obtain MFs protected by the mirror symmetry. 

Among the two possible mirror symmetries A4 xy in 
Eq. only M% y satisfies Eq. ©. Indeed, from 

A4 xy = —1 = A 2 , one can show that A = ±i and 
(M%Y = X 2 M^, which yields CM±p = ±X 2 M±* 



'xy ■ 



These are consistent with Eq. ([5]) only for Ai xy , which 
means that when the gap function is even under the mir- 
ror reflection AA xy A(k)A4 xy = A(fc), the subsector with 
the eigenvalue A supports its own particle-hole symme- 
try [Fig. [Ha)]. However, when the gap function is odd 
under the mirror reflection M xy A(k)M t xy = -A(fc), 
the particle-hole symmetry of the original Hamiltonian 
merely interchanges two different subsectors, and thus 
each subsector does not have particle-hole symmetry 

[Fig. mb)]. 

In the latter case, we find that the subsectors belong 
to class A of the topological classification 2l| since there 
is no particle-hole symmetry. That is, they are topologi- 
cally the same as quantum Hall states, which also belong 
to class A. Therefore, even when there are topologically 
protected states ensured by the mirror Chern number, 
these states reduce to Dirac fermions as is the case for 
quantum Hall states. 

In contrast, MFs can be realized in the former case: 
each subsector belongs to class D of the topological clas- 
sification 21[ because of its own particle-hole symmetry. 
Consequently, as a class D topological phase, the mir- 
ror subsector can host MFs. The topological invariant 
for class D is for three dimensions, Z for two dimen- 
sions, and Zi for one dimension. Therefore, in addition 
to the mirror Chern number u(X), we can define the one- 
dimensional (ID) Z2 invariants ^(A) and v y (X) as Q 



v x {X) = - 

7T 



dk y A y (ir, k y ) mod 2, 



(7) 



which we refer to as mirror Zi invariants (y y {X) is de- 
fined in a similar manner). The ID mirror Z2 invariants 
are well defined even for a system without time-reversal 
invariance, and thus they arc very different from the Zi 
invariants for topological (crystalline) insulators. 

The topological invariants ^(A), ^(A), and ^(A) char- 
acterize the MFs in the mirror subsector; the mirror 
Chcrn number ensures the existence of chiral edge MFs 
in the mirror subsector. U sing a similar argument to 
that given in Refs. [13 and [31(, we can also show that 
if the mirror Chern number is odd, then the Majorana 
zero-energy states exist on a vortex. Furthermore, the ID 
mirror Zi invariants guarantee the existence of Majorana 
zero-energy bound states localized on a dislocation, if the 
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Majorana fermion 
(a) / , 1. 




FIG. 1: Two possible realizations of particle-hole symmetry. 



Burgers vector B characterizing the dislocation satisfies 
B G x = 1 mod 2 [32]. Here, G A = ±(v x (X)b x +i/y(\)b y ) 
with b x and b y being the reciprocal lattice vectors in the 
x- and y-directions. All these MFs are topologically sta- 
ble so long as the mirror symmetry is preserved macro- 
scopically beyond the scale of the coherence length. 

A generalization of these results to three-dimensional 
SCs is fairly straightforward. In three dimensions, 
Eq. J2]) is replaced with M xy H(k x , k y , k z )M xy = 
H(k x ,k y , — fc 2 ). If we then consider the mirror invariant 
fczfcy-planes with k z = or n, the same argument as per 
the two-dimensional case applies, and again symmetry- 
protected MFs can be realized when the gap function is 
even under the mirror reflection. We can also general- 
ize our arguments to other point group symmetries. For 
instance, in a manner similar to that used for topolog- 
ical crystalline insulators, we can introduce topological 
invariants using discrete rotation symmetries, but to sta- 
bilize the MFs, additional requirements similar to those 
in Eq. (J6|) arc needed. 

To illustrate the general arguments above, we ap- 
ply our results to the two-dimensional spin-triplet SC 
Sr2Ru04 [HI, Hi)]. Sr2RuC>4 has a tetragonal structure 
with the crystal point group symmetry D±h and, in par- 
ticular, is invariant under a mirror reflection with respect 
to the iry-plane. We consider the superconducting state 
of Sr2Ru04 under magnetic fields in the z-direction. In 
this case, the crystal point group symmetry reduces to 
Cih , but the system is still invariant under the mirror re- 
flection A4 xy . According to the representation theory of 
C4/J, the possible spin-triplet gap functions are classified 
as in Table |U 

For zero or weak magnetic fields, a number of experi- 
ments support time-reversal breaking chiral p x + ip y su- 
perconductivity in Sr2Ru04, which belongs to the E u 
representation in Table HI (33. 34 1. The d- vector is aligned 
along the ^-direction, and the gap function, A(fc) = 



id(k) ■ cra y , is odd under the mirror reflection. For mag- 
netic fields greater than 20 mT, Knight shift measure- 
ments suggest that the d-vector is parallel to the xy- 
plane, consistent with the A u or B u representation in 
Table U (35[. Thus, the gap function is helical and even 
under the mirror reflection. For consistency with both 
experimental results, there should be a phase transition 
at a critical magnetic field H zc at which the d-vector 
rotates from the z-direction to the xy-plane, the chiral 



TABLE I: Possible ei-vector states of Sr2RuC>4 under mag- 
netic fields in the z-direction. The representation F, d-vector, 
parity of the gap function A = id ■ cra y under the mirror re- 
flection M X y, topological class of the subsector with a definite 
eigenvalue of M xy , and mirror topological numbers v, v x , and 
v y are summarized. Note that the mirror Z-2 numbers are not 
defined in the E u representation. 



r 


d-vector 


M xy 


subsector 


MA),^(A), 


u y (X)] 


A u 


x sin k x + y sin k y 
x sin k y — y sin k x 


even 


class D 


L±1,1,1Ja 


=±i 


B u 


x sin k x — V sin k y 
x sin k y + y sin k x 


even 


class D 


[=fM,i]a 


=±i 


E u 


z(sin k x + i sin k y ) 


odd 


class A 


[1,- Ha 


=±i 



superconducting state becomes helical, and the mirror 
parity of the gap function changes from odd to even. 

It is known that the chiral p x + ip y superconducting 
state in the low-field phase hosts topological supercon- 
ductivity [H, [37|. Taking the spin degrees of freedom 
into account, the Chern number of this state is evaluated 
as i^ch = 2, and thus a pair of topologic ally protected 
chiral fermions exist on the boundary However, 
our argument implies that it is difficult to detect their 
Majorana character. Since the gap function is odd under 
the mirror reflection, this state can be separated into a 
pair of subsectors belonging to class A. In each subsec- 
tor, the paired MFs, 71 and 72, are recast into a single 
Dirac fermion i\) = 71 + 472, and thus their dynamics are 
reduced to those of the Dirac fermion. 



In contrast, the helical superconducting state in the 
high-field phase supports MFs. The ordinary Chern num- 
ber z^ch is found to be zero in this phase, and so this phase 
is not topological in terms of the topological periodic ta- 
ble [2l|. Nevertheless, since the gap function is even un- 
der the mirror reflection, the mirror Chern number and 
the ID mirror Z2 invariants can be introduced. Using a 
standard adiabatic deformation method for the Hamilto- 
nian, the mirror topological invariants can be calculated, 
and we find that they are nonzero. The obtained results 
are summarized in Table HI From the bulk-boundary cor- 
respondence, one can conclude that MFs protected by the 
mirror symmetry exist. Since the mirror symmetry is es- 
sential for hosting the topological phase, the high-field 
phase realizes a topological crystalline SC[44|. 

To illustrate these results, we examine the quasiparti- 
cle states of Si'2Ru04 using the following model Hamil- 
tonian on the square lattice [HI . The conduction bands 
of Sr 2 Ru04 consist of three 4d-£2 9 Ru orbitals, d xz , d yz , 
and d xy , which we label as I = 1,2, and 3, respectively. 
The superconducting state of Sr2Ru04 is described by 



4 



ClassA(M,, = i) 




index of states 



Class D(M„= i) 



(b) - 1 




index of states 



dllz Low filed phase 



d-Lz High field phase 



->• Hz 



FIG. 2: Topological gapless states in the low- and high-field 
phases of S^RuCU. (a) and (b) Edge states and (c) and (d) 
bound states of the edge dislocation. In (a) and (b), the red 
(blue) lines represent edge states localized at x = [x = L). 



Results are shown for the M x 



i sector. The d-vector is 



chosen as d(k) = z(sink x + isink y ) in the low-field phase 
and as d(k) — x sin k y — y sin k x in the high-field phase. We 
take ti=t 3 = 0.5, t 3 = 0.2, U = 0.1, it = -0.2, ti' = -0.2, 



HbH z = 0.1, A = 0.3, and A 



0.' 



the BdG Hamiltonian H = %kin + %so + %pair, in which 



£fci 9k 
H kin - ^(4 sl ,4 s2 ,4 s3 ) g k e k2 

ks \ £fc 3 




(8) 



kiss' 



where e^i = —2^ cos fcj, — /i, = — 2ij cos k x — fj,, and 
£fc3 = —2t 2 {cosk x + cosfc.y) — 4i 3 cos fcj, cos fc y — fi' are 
the kinetic terms of orbital 1, 2, and 3, respectively, 
g k = — 4i4 sin fc^ sin fc y is the hybridization of orbital 1 
and 2, and A'(fc) = iA l d(k) ■ <ra y is the gap func- 
tion with d(k) as defined in Tabic |TJ Here, <t" is the 
Pauli matrix, and ei mn is the completely antisymmetric 
tensor. In the (c ks i,c ks2 ,c ks3 ) basis, M xy is given by 
■M X y = diag(— icr 2 , — i<r z , io z ). For definitcness, we as- 
sume d(k) = x sin k y — y sin k x in the high-field phase, 
but a qualitatively similar result is obtained as long as 
the d-vector is parallel to the ccy-plane. 

To examine the edge states of Sr 2 Ru04, we numeri- 
cally diagonalize the BdG Hamiltonian Eq. ([5]) in the 
coordinate space with open boundaries at x = and 
x = L = 30. Using the mirror symmetry, we separate 
the quasiparticle spectra into two subsectors with differ- 
ent eigenvalues of M. xy . Figures (2^ a) and (b) illustrate 
the quasiparticle spectra in the M. xy — i sector. A sim- 
ilar result is obtained in the A4 xy = —i sector. Gapless 
edge states for both the low-field phase [Fig. [3] (a)] and 



the high-field phase [Fig. [2](b)] can be seen clearly, which 
is consistent with the nonzero mirror Chern numbers in 
both phases. Note, however, that there is an essential dif- 
ference between the two phases. While the quasiparticle 
spectrum in the high-field phase is particle-hole symmet- 
ric as E(k y ) = —E(—k y ), that in the low-field phase is 
not. This indicates clearly that only the high-field phase 
supports edge MFs. Note also that the particle-hole sym- 
metry in the low-field phase is recovered if it is considered 
in conjunction with the Ai xy = —i sector. 

The bound states of the edge dislocations are shown 
in Figs. ^c) and (d). The energy spectra were calcu- 
lated on 30 x 30 unit cell system with periodic boundary 
conditions in the x- and y-dircctions. Two edge dislo- 
cations with the Burgers vector B = ±e x , separated by 
half the length of the system size, are considered. In the 
high-field phase, we have G x = j^{b x + b y ) from Ta- 
bled! so G x ■ B = 1, which predicts the existence of a 
zero mode in the dislocation. The high-field phase can 
be seen to support zero energy states, which are local- 
ized on the edge dislocations and are two-fold degenerate 
because there are two dislocations in the system. There- 
fore, there is a single Majorana zero mode localized on 
each edge dislocation. 

In the same way as ordinary MFs, these symmetry- 
protected MFs can be identified by tunneling spec- 
troscopy. In particular, the detection of the zero mode 
in the dislocation gives a strong signature of the d- vector 
rotation in the high-field phase; if the d-vector rotation 
does not occur, the bound state in the dislocation must 
have a gap of 0{hbH z ). Thus, the observation of the 
zero energy state is a distinct signal of the d-vector ro- 
tation. The change in the edge-mode spectra can also 
provide details of the tunneling conductance. In addi- 
tion, in each sector of a definite eigenvalue of A4 xy , the 
symmetry-protected MFs behave like ordinary MFs and 
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retain the unique features specific to MFs [46 

The symmetry-protected MFs discussed here are appli- 
cable to many unconventional SC; the same arguments 
apply equally to Cu^Bi^es SBO] and UPt 3 |H,|6^ 
as both these crystals have a mirror plane. Details about 
these systems will be reported elsewhere. 
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